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We present the measurement of non-photonic electron production at high transverse momentum 
(pr > 2.5 GeV/c) in p+p collisions at ^/s = 200 GeV using data recorded during 2005 and 2008 by 
the STAR experiment at the Relativistic Heavy Ion Collider (RHIC). The measured cross-sections 
from the two runs are consistent with each other despite a large difference in photonic background 
levels due to different detector configurations. We compare the measured non-photonic electron 
cross-sections with previously published RHIC data and pQCD calculations. Using the relative con- 
tributions of B and D mesons to non-photonic electrons, we determine the integrated cross sections 
of electrons ( e ) at 3 GeV/c < pr < 10 GeV/c from bottom and charm meson decays to be 
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(B ^ e) +<f^ p ^ e ? \ ye=0 = 4.0±0.5(stat.)±l.l(syst.) nb and 
nb, respectively. 

PACS numbers: 13.20.Fc, 13.20.Hc, 25.75.Cj 



,=o = 6.2±0.7(stat.)±1.5(syst. 



I. INTRODUCTION 

Heavy quark production in high-energy hadronic colli- 
sions has been a focus of interest for years. It is one of 
the few instances in which experimental measurements 
can be compared with QCD predictions over nearly the 
entire kincmatical range nH3f. Due to the large masses 
of charm and bottom quarks, they are produced almost 
exclusively during the initial high-Q parton-parton inter- 
actions and thus can be described by perturbative QCD 
calculations. 

Measurement of heavy flavor production in elementary 
collisions represents a crucial test of the validity of the 
current theoretical framework and its phenomenological 
inputs. It is also mandatory as a baseline for the inter- 
pretation of heavy flavor production in nucleus-nucleus 
collisions In these heavy-ion collisions one investi- 
gates the properties of the quark-gluon plasma (QGP), 
which is created at sufficiently high center-of-mass ener- 
gies. Many effects on heavy flavor production in heavy- 
ion collisions have been observed but are quantitatively 
not yet fully understood Q. Of particular interest are 
effects which modify the transverse momentum spectra 
of heavy flavor hadrons, including energy loss in the QGP 



("jet quenching" 
as elliptic flow LL 




as well as collective effects such 
In addition, J/iJj might be regen- 
erated in a dense plasma from the initial open charm yield 
fl2| , making precise measurements of the transverse mo- 
mentum spectra in elementary p+p collisions imperative. 

Open heavy-flavor production in p+p, d+A, and A+A 
collisions at ^/snn = 200 GeV has been studied at the 
Relativistic Heavy Ion Collider (RHIC) using a variety 
of final-state observables Q. The STAR collaboration 
measured charm mesons directly through their hadronic 
decay channels fl3l - ll5l |. Due to the lack of precise 
secondary vertex tracking and trigger capabilities these 
measurements are restricted to low momenta {pt < 3 
GeV/c). Both STAR QJ, d and PHENIX 0,13 also 
measured heavy flavor production through semileptonic 
decays of charm and bottom mesons (D,B -+ I vg X\ 
While the measured decay leptons provide only limited 
information on the original kinematics of the heavy flavor 
parton, these measurements are facilitated by fast online 
triggers and extend the kinematic range to high px- 

In this paper, we report STAR results on non-photonic 
electron production at midrapidity in p+p collisions at 
yfs = 200 GcV using data recorded during year 2005 
(Run2005) and year 2008 (Run2008) with a total in- 
tegrated luminosity of 2.8 pb _1 and 2.6 pb _1 , respec- 
tively. The present results are consistent with the Next- 
to-Leading Logarithm (FONLL) calculation within its 
theoretical uncertainties. Utilizing the measured rela- 
tive contributions of B and D mesons to non-photonic 



electrons which were obtained from a study of electron- 
hadron correlations (e-h) [l9j , we determine the invariant 
cross section of electrons from bottom and charm meson 
decays separately at pr > 3.0 GeV/c. 

The article is organized as follows. In Sec. |H] we de- 
scribe the STAR detectors and triggers relevant to this 
analysis. Sec. IIIII describes the data analysis in detail, 
and in Sec. HVI we present and discuss the results. Sec. W\ 
provides conclusions. 



II. EXPERIMENT 

A. Detectors 

STAR is a large acceptance, multi-purpose experiment 
composed of several individual detector subsystems with 
tracking inside a large solenoidal magnet generating a 
uniform field of 0.5 T [2(|. The detector subsystems rel- 
evant for the present analysis are briefly described in the 
following. 



1. Time Projection Chamber 

The Time Projection Chamber (TPC) HH is the main 
charged particle tracking device in STAR. The TPC cov- 
ers ±1.0 units in pscudorapidity (rj) for tracks crossing all 
layers of pads, and the full azimuth. Particle momentum 
is determined from track curvature in the solenoidal field. 
In this analysis, TPC tracks are used for momentum de- 
termination, electron-hadron separation (using specific 
ionization dE/dx), to reconstruct the interaction vertex, 
and to project to the calorimeter for further hadron re- 
jection. 



Barrel Electromagnetic Calorimeter and 
Barrel Shower Maximum Detector 



The Barrel Electromagnetic Calorimeter (BEMC) 
measures the energy deposited by photons and electrons 
and provides a trigger signal. It is located inside the 
magnet coil outside the TPC, covering | ^7 1 < 1.0 and 
27r in azimuth, matching the TPC acceptance. The 
BEMC is a lead-scintillator sampling electromagnetic 
calorimeter with a nominal energy resolution of 5E /E ~ 
14%/^E/lGeV © 1.5% [12]. The full calorimeter is 
segmented into 4800 projective towers. A tower covers 
0.05 rad in (j> and 0.05 units in r\. Each tower consists of 
a stack of 20 layers of lead and 21 layers of scintillator 
with an active depth of 23.5 cm. The first two scintillator 
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layers are read out separately providing the calorimeter 
preshower signal, which is not used in this analysis. A 
Shower Maximum Detector (BSMD) is positioned behind 
the fifth scintillator layer. The BSMD is a double layer 
wire proportional counter with strip readout. The two 
layers of the BSMD, each containing 18000 strips, pro- 
vide precise spacial resolution in <f> and r\ and improve 
the electron-hadron separation. The BEMC also pro- 
vides a high-energy trigger based on the highest energy 
measured by a single tower in order to enrich the event 
samples with high-p^ electromagnetic energy deposition. 



require the radial location of the first TPC point of re- 
constructed tracks to be less than 70 cm (see Sec. IIII Bp 
in the Run2008 analysis; furthermore, conversion elec- 
trons originating from TPC gas have low probability to 
be reconstructed by the TPC tracking due to the short 
track length. While the Run2008 simulation describes 
the material distribution very well, the material budget 
for the support structure and electronics related to the 
silicon detectors is not reliably described in the Run2005 
simulation [27}. This, however, has little effect on this 
analysis, as explained in Sec. IIIIDl 



3. Trigger Detectors 

The Beam-Beam Counters (BBC) [23| are two identi- 
cal counters located on each side of the interaction region 
covering the full azimuth and 2.1 < \q\ <5.0. Each de- 
tector consists of sets of small and large hexagonal scin- 
tillator tiles grouped into a ring and mounted around the 
beam pipe at a distance of 3.7 m from the interaction 
point. In both Run2008 and Run2005, the BBC served 
as a minimum-bias trigger to record the integrated lumi- 
nosity by requiring a coincidence of signals in at least one 
of the small tiles (3.3 < |ry| <5.0) on each side of the inter- 
action region. The cross-section sampled with the BBC 
trigger is 26.1 ± 0.2(stat.) ± 1.8(syst.) mb [H| for p+p 
collisions. The timing signal recorded by the two BBC 
counters can be used to reconstruct the collision vertex 
along the beam direction with an accuracy of ~ 40 cm. 

The data in d+Au collisions recorded during year 2008 
is used as a crosscheck in this analysis (see Sec IIII El) . 
During this run, a pair of Vertex Position Detectors 
(VPD) [25[ was also used to select events. Each VPD 
consists of 19 lead converters plus plastic scintillators 
with photomultiplier-tube readout that are positioned 
very close to the beam pipe on each side of STAR. Each 
VPD is approximately 5.7 m from the interaction point 
and covers the pseudorapidity interval 4.24 < \r/\ < 5.1. 
The VPD trigger condition is similar to that of the BBC 
trigger except that the VPD has much better timing res- 
olution, enabling the selected events to be constrained 
to a smaller range (~ ±30 cm in d+Au run) around the 
interaction point. 



B. Material Thickness in front of the TPC 

Table [T] shows a rough estimate of material thickness 
between the interaction point and the inner field cage 
(IFC) of the TPC during Run2008 in the region relevant 
to the analysis. The amount of material is mostly from 
the beam pipe (BP), the IFC, air, and a wrap around 
the beam pipe. In Run2005 , the amount of material is 
estimated to be ~ 10 times larger in front of the TPC 
inner field cage [26| and is dominated by the silicon de- 
tectors which were removed before Run2008. The con- 
tribution from the TPC gas is not significant because we 



TABLE I: Estimates of material thickness of the beam pipe, 
the wrap around the beam pipe, the TPC inner field cage 
and air between the beam pipe and the inner field cage in 
Run2008. 



source 


thickness in radiation lengths 


beam pipe 


0.29 % 


beam pipe wrap 


~ 0.14 % 


aii- 


~ 0.17 % 


inner field cage 


~ 0.45 % 



C. Triggers and Datasets 

The data reported in this paper were recorded during 
Run2005 and Run2008 at y/s= 200 GeV. All events used 
in this analysis are required to satisfy a BEMC trigger 
and a BBC minimum-bias trigger. In addition, event 
samples using a VPD trigger in the 2008 d+Au colli- 
sions are used for systematic cross-checks as described in 

secnm 

To enrich the data sample with high-px electromag- 
netic energy deposition, the BEMC trigger requires the 
energy deposition in at least one tower to exceed a pre- 
set threshold (high-tower). Most of the energy from 
an electron or a photon will be deposited into a sin- 
gle tower since the tower size exceeds the radius of a 
typical electromagnetic shower. The Run2008 datasets 
used here were recorded using three high-tower triggers 
with different thresholds, corresponding to a sampled lu- 
minosity of ^2.6 pb _1 . Expressed in terms of trans- 
verse energy (Et), the thresholds were approximately 
2.6 GeV, 3.6 GeV and 4.3 GeV. The Run2005 datasets 
used here are from two high-tower triggers with Et 
thresholds of 2.6 GeV (HT1) and 3.5 GeV (HT2), cor- 
responding to a sampled luminosity of ^2.8 pb _1 . In 
the Run2008 analysis, datasets from different high-tower 
triggers are treated together after being combined, with 
double counting avoided by removing duplicates in the 
corresponding high-tower ADC spectra. Trigger efficien- 
cies and prescale factors imposed by the data acquisition 
system are taken into account during the combination. 
In the analysis of the Run2005 data, HT1 and HT2 data 
are treated separately 
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In Run2005 the integrated luminosity was monitored 
using the BBC minimum-bias trigger, while in Run2008, 
because of the large beam related background due to high 
luminosity, a high threshold high-tower trigger seeing a 
total cross section of 1.49 fjb, was used as luminosity 
monitor. 



III. ANALYSIS 
A. Photonic Background Removal 

The main background in this analysis is the substantial 
flux of photonic electrons from photon conversion in the 
detector material and Dalitz decay of ir° and r\ mesons. 
These contributions need to be subtracted in order to 
extract the non-photonic electron yield, which is dom- 
inated by electrons from semileptonic decays of heavy 
flavor mesons. 

There are two distinct methods for evaluating contri- 
butions from photonic electrons. In the cocktail method, 
the estimated or measured invariant cross-sections are 
used to calculate contributions from various sources 
(mostly 7T°, r\ mesons), and to derive from those the pho- 
tonic electron distributions. Given sufficient knowledge 
of the production yield of those mesons, this method al- 
lows one to determine directly the contributions from 
Dalitz decays. With this method, a detailed understand- 
ing of the material distribution in the detector is re- 
quired in order to evaluate the contribution from pho- 
ton conversion. Another method, used in this analysis, 
is less dependent on the exact knowledge of the amount 
of material. This method reconstructs the photonic elec- 
trons through the specific feature that photonic electron- 
positron pairs have very small invariant mass. Not all 
photonic electrons can be identified this way since one of 
the electrons may fall outside of the detector acceptance, 
or has a very low momentum, in which cases both elec- 
trons in the pair are not reconstructed. This inefficiency 
must be estimated through simulation. 

Electron pairs are formed by combining an electron 
with pt > 2.5 GcV/c, which we refer to as a primary 
electron, with all other electrons (partners) reconstructed 
in the same event, with opposite charge sign (unlikc- 
sign) or same charge sign (like-sign). The upper two 
panels of Fig. [T] show the invariant mass spectra for pri- 
mary electrons with 2.5 GeV/c < px < 3.0 GeV/c (left) 
and 8.0 GcV/c < p T < 10.0 GcV/c. The unlike-sign 
spectrum includes pairs originating from photon conver- 
sion and Dalitz decay, as well as combinatorial back- 
ground. The latter can be estimated using the like-sign 
pair spectrum. The photonic electron spectrum is ob- 
tained by subtracting like-sign from unlike-sign spectrum 
(unlike-minus-like) . The broad shoulders extending to- 
ward higher masses in the spectra are caused by finite 
tracking resolution, which leads to a larger reconstructed 
opening angle when the reconstructed track helices of two 
conversion electrons intersect each other in the transverse 
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FIG. 1: (Color online) The upper two panels show the electron 
pair invariant mass distributions for electrons at 2.5 GeV/c < 
p T < 3.0 GeV/c (a) and at 8.0 GeV/c < p T < 10.0 GeV/c 
(b). Solid and dashed lines represent unlike-sign and like-sign 
pairs, respectively. Closed circles represent the difference of 
unlike and like. The lower two panels show the simulated 
invariant mass spectra with electrons at 2.5 GeV/c < pr < 
3.0 GeV/c (c) and at 8.0 GeV/c < p T < 10.0 GeV/c (d). 
Solid and dashed lines represent results from 7 conversions 
and 7T° Dalitz decay. 

plane. The overall width of the mass spectra depends on 
the primary electron pt, but most photonic pairs are 
contained in range of m ee < 0.24 GeV/c 2 . The lower 
two panels of Fig. [T] show the simulated invariant mass 
spectra of the two dominant sources of photonic elec- 
trons, 7T° Dalitz and 7 conversions, in the same two pt 
regions, which are similar in shape due to similar decay 
kinematics. The electron spectrum obtained from the 
unlike-minus-like method is from pure photonic electrons 
because the combinatorial background is accurately de- 
scribed by the like-sign pair spectrum according to the 
simulation and the simulated mass spectra are in quali- 
tative agreement with the data. This is also proved by 
the fact that the distribution of the normalized ioniza- 
tion energy loss (see Sec. IIII B[) can be well described by 
a Gaussian function expected from electrons as shown in 
Fig.|H 

We calculate the yield of non-photonic electrons ac- 
cording to 

N(npe) = N(inc)- e purlty - N(pho)/e pho , 

where N(npe) is the non-photonic electron yield, N(inc) 
is the inclusive electron yield, N(pho) is the photonic 
electron yield, e p h is the photonic electron reconstruc- 
tion efficiency defined as the fraction of the photonic elec- 
trons identified through invariant mass reconstruction, 
and €p Ur ity is the purity reflecting hadron contamination 
in the inclusive electron sample. 
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FIG. 2: (a) The distribution of the minimum distance be- 
tween an electron track projection point at the BEMC and 
all BEMC clusters along the beam direction from unlike-sign 
electron candidate pairs (solid line), like-sign electron can- 
didate pairs (dashed line) and unlike-minus-like (closed cir- 
cles), (b) p/Eo distribution from unlike-sign electron can- 
didate pairs (solid line), like-sign electron candidate pairs 
(dashed line) and unlike-minus-like (closed circles). 



Electrons from open heavy flavor decays dominate non- 
photonic electrons. The contribution from semi-leptonic 
decay of kaons is negligible for pr > 2.5 GcV/c fl7| |. 
Electrons from vector mesons (p,u),<f>, J/i/j, T) decays 
and Drell-Yan processes are subtracted from the mea- 
surement (see Sec. IIII Gl for details). 



B. Electron Reconstruction and Identification 
Efficiency 

In the analysis of the Run2008 data, we select only 
tracks with pr > 2.5 GeV/c and |?7| < 0.5. The event 
vertex position along the beam-axis (V z ) is required to 
be close to the center of the TPC, i.e. \V Z \ < 30 cm. To 
avoid track reconstruction artifacts, such as track split- 
ting, the selected tracks are required to have at least 52% 
of the maximum number of TPC points allowed in the 
TPC, a minimum of 20 TPC points and a distance-of- 
closest-approach (DCA) to the collision vertex less than 
1.5 cm. For hadron rejection we apply a cut of na e > — 1 
on the normalized ionization energy loss in the TPC [28j . 
which is defined as 



ncr e 



log{{dE/dx)/B e )/a e 



where B e is the expected mean dE/dx of an electron 
calculated from the Bichsel function [29} and a e is the 
TPC resolution of \og{{dE/dx)/B e ). 

We reconstruct clusters in the BEMC and the BSMD 
by grouping adjacent hits that are likely to have origi- 
nated from the same incident particle [27|. The selected 
tracks are extrapolated to the BEMC and the BSMD 
where they are associated with the closest clusters. The 



association windows for electrons are determined by mea- 
suring the distance between the track projection point 
at the BEMC (BSMD) and the closest BEMC (BSMD) 
cluster using photonic electrons from the unlike-minus- 
like pairs. Figure [2] (a) shows the distribution of this 
distance at the BEMC along the beam direction for elec- 
trons from unlike-sign, like-sign and unlike-minus-like 
pairs requiring m e + e - < 0.24 GeV/c 2 , a maximum 1.0 
cm DCA between two helical-shaped electron tracks and 
a 3.0 keV/cm < dE/dx < 5.0 keV/cm cut on ionization 
energy loss for partner tracks. Most electrons are inside 
a window of ±20 cm around the track projection point. 
The window in the azimuthal plane is determined to be 
±0.2 radian. Figurc[2](b) shows the distribution oip/Eo 
for electrons from unlike-sign, like-sign and unlike-minus- 
like pairs, where Eq is the energy of the most energetic 
tower in a BEMC cluster. The distribution is peaked 
around one due to the small mass of the electron and the 
fact that most electron energy is deposited into one tower. 
We apply a cut of p/E < 2.0 to further reduce hadron 
contamination. Cuts on the association with BSMD clus- 
ters are kept loose to maintain high efficiency. Each track 
is required to have more than one associated BSMD strip 
in both <f> and r\ planes. 

The efficiencies of electron identification cuts are esti- 
mated directly from the data using pure photonic elec- 
trons obtained from the unlike-minus-like pairs requiring 
m e + e - < 0.24 GeV/c 2 , a maximum 1.0 cm DCA between 
two helical-shaped electron tracks and a 3.0 keV/cm 
< dE/dx < 5.0 keV/cm cut on ionization energy loss 
for partner tracks. A cut of pr > 0.3 GeV/c for partners 
is also applied, selecting a region where the simulation 
does a good job of describing the data. The efficiency for 
one specific cut is then calculated as the ratio of electron 
yield before the cut to that after the cut, while all the 
other electron identification cuts are applied. To avoid 
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FIG. 3: (Color online) Efficiencies of the cuts on number of 
TPC points (open circles), na e (open triangles) and BEMC 
(open squares) in (a) the Run2008 and (b) the Run2005 anal- 
yses. The closed circles represent the total efficiency which is 
the product of all individual ones. 
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possible correlation among different cuts, efficiencies for 
all BEMC and BSMD cuts are calculated together. Fig- 
ure |3] (a) shows the breakdown of the electron identifica- 
tion efficiency as a function of pr ■ The drop in the low 
Pt region comes mainly from BSMD inefficiency, while 
the drop in the high pr region is caused by the p/Eq cut. 
The uncertainties in the figure are purely statistical and 
are included as part of the systematic uncertainties for 
the cross section calculation. 

To maintain high track quality and suppress photonic 
electrons from conversion in the TPC gas, we require the 
radial location of the first TPC point to be less than 70 
cm. This cut causes an inefficiency of 12.0 ± 2.0% for 
non-photonic electrons according to the estimates from 
both simulation and data. 

Through embedding simulated electrons into high- 
tower trigger events and processing them through the 
same software used for data production, single elec- 
tron reconstruction efficiency in the TPC is found to be 

0. 84±0.04 with little dependence on momentum for px > 
2.0 GcV/c. 

The analysis of the Run2005 data is slightly differ- 
ent from that of the Run2008 described above. Only 
half of the BEMC (0 < rj < 1.0) was instrumented in 
2005. Due to the presence of the silicon detectors, and 
their significant material budget, photonic backgrounds 
were substantially higher. We select only tracks with 
< rj < 0.5 from -30 cm < V z < 20 cm in order to 
avoid the supporting cone for the silicon detectors in the 
fiducial volume while keeping track quality cuts identi- 
cal to those in the Run2008 analysis. However, we apply 
a tighter cut on the normalized ionization energy loss, 

1. e. —0.7 < na e < 3.0, to improve hadron rejection. 
BEMC clusters are grouped with geometrically overlap- 
ping BSMD clusters to improve position resolution and 
electron hadron discrimination through shower profile. 
The clustering algorithm is also modified to increase the 
efficiency of differentiating two overlapping BSMD clus- 
ters by lowering the energy threshold of the second cluster 
[30l ] . The minimum angle between track projection point 
at the BEMC and all BEMC clusters is required to be less 
than 0.05 radian. We also require each track to have more 
than one associated BSMD strip in both <p and r\ planes, 
and a tightened p/E cut of 0.3 < p/E < 1.5, where E is 
the energy of the associated BEMC cluster. The efficien- 
cies for the electron identification cuts are estimated by 
embeddin g si mulated single electrons into minimum-bias 
PYTHIA [31j events. Figure [3] (b) shows the breakdown 
of electron identification efficiency as a function of pt 
in the Run2005 analysis. There is no drop at high pt 
as in the Run2008 result because the energy of a whole 
BEMC cluster, instead of the highest tower, is used for 
the p/E cut. No cut on the first TPC point is applied in 
this analysis. To avoid the TPC tracking resolution effect 
that causes the broad shoulder extending toward higher 
masses in the invariant mass spectrum of the Run2008 
analysis, we utilize a 2-D invariant mass by ignoring the 
opening angle in the 4> plane when reconstructing the 
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FIG. 4: (Color online) Left two panels are na e distributions 
in the Run2008 analysis for unlike-sign (dot-dashed line) , like- 
sign (dotted line) and unlike-minus-like (closed circles) pairs 
together with a Gaussian fit (solid lines) at (a) 2.5 GeV/c < 
p T < 3.0 GeV/c and (b) 8.0 GeV/c < p T < 10.0 GeV/c 
after applying all the electron identification cuts except the 
na e cut. Right panel (c) shows the mean and width of the 
Gaussian fitting functions for pure photonic electron (unlike- 
minus-like) n<j e distribution as shown in left panels for each 
Pt bin. See text for details. 



e + e~ invariant mass [301 ] - We require —3 < na e < 3 
for partner tracks, 2-D m e + e - < 0.1 GeV/c 2 for pairs, 
a maximum 0.1/0.05 radian for the opening angle in the 
4>/9 plane, and a maximum 1.0 cm DC A between two 
electron helices. A cut of pr > 0.3 GeV/c for partners is 
also applied so that the simulation can describe the data 
well. 

By following independent analysis procedures from two 
RHIC runs where the amount of material for photonic 
background is significantly different, we will be able to 
validate our approach for measuring non-photonic elec- 
tron production. 



C. Purity Estimation 

After applying all electron identification cuts, the in- 
clusive sample of primary electrons is still contaminated 
with hadrons. To estimate the purity of electrons in the 
inclusive sample, we perform a constrained fit on the 
charged track na e distributions in different pt regions 
with three Gaussian functions representing the expected 
distributions of tp* 1 , K +p and e . The purity is es- 
timated from the fit. 

Ideally the electron na e will follow the standard nor- 
mal distribution. The actual distribution can be slightly 
different due to various effects in data calibrations. We 
can, however, determine its shape in different px re- 
gions directly from data using photonic electrons from 
the unlike-minus-like pairs. The left panel of Fig. |4] shows 
the ncr e distribution for tracks with (a) 2.5 GeV/c < 
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FIG. 5: (Color online) na e distribution for inclusive electrons 
(closed circles) and fits from different components at (a) 2.5 
GeV/c < p T < 3.0 GeV/c in the Run2008 analysis, (b) 2.5 
GeV/c < pr < 3.5 GeV/c in the Run2005 analysis and (c) 
8.0 GeV/c <p T < 10. GeV/c in the Run2008 analysis after 
applying all electron identification cuts except the na e cut. 
Different curves represent K + p (dotted line), -k^ 1 (dot- 
dashed line), electrons (dashed line) and the overall fit (solid 
lines) 



p T < 3.0 GeV/c and (b) 8.0 GcV/c < p T < 10.0 GeV/c 
from unlike-sign, like-sign pairs as well as for photonic 
electrons from the unlike- minus- like pairs. Here all elec- 
tron identification cuts, except the na e cut, are applied. 
The na e of photonic electrons are well fitted with Gaus- 
sian functions. Figure U (c) shows the mean and width 
of the Gaussian fit as a function of electron px, which, 
as discussed above, differ slightly from the ideal values. 
The solid lines in the figure are fits to the data using 
a second order polynomial function. The dotted lines 
are also second order polynomial fits to the data except 
that the data points are moved up and down simulta- 
neously by one standard deviation. The region between 
the dotted lines represents a conservative estimate of the 
fit uncertainty since we assume that the points are fully 
correlated. The mean, width and their corresponding un- 
certainties from the fits are used to define the shape of 
electron na e distribution in the following 3- Gaussian fit. 
The na E of tt^ 1 and + p ± are also expected to follow 
Gaussian distributions |28j. Ideally their width is one 
and their means can be calculated through the Bichsel 
function [2{|. These ideal values are used as the initial 
values of the fit parameters in the following 3-Gaussian 
fit. 

Figure [5] shows the constrained 3-Gaussian fits to the 
na e distributions of inclusive electron candidates with 
2.5 GeV/c < p T < 3.0 GeV/c in the Run2008 analy- 
sis (upper-left), 2.5 GcV/c < pt < 3.5 GcV/c in the 
Run2005 analysis (lower-left) and 8.0 GeV/c < pt < 
10.0 GeV/c in the Run2008 analysis (right). Here we 
leave the na e cut open. The dotted, dot-dashed and 



FIG. 6: Purity of the inclusive electron sample as a function 
of p T in data from (a) Run2008 and (b) Run2005. The 2008 
result is from combined datasets of all different high-tower 
triggers. The 2005 results for the two different high-tower 
triggers, i.e. HT1 (closed circles) and HT2 (open circles), are 
plotted separately. 



dashed lines represent, respectively, the fits for K ± +p , 
7T , and e . Compared to the Run2008 analysis, the 
electron component in the Run2005 analysis at similar 
Pt is more prominent due to the larger conversion elec- 
tron yield. The solid lines are the overall fits to the spec- 
tra. The purity is calculated as the ratio of the integral 
of the electron fit function to that of the overall fit func- 
tion above the na e cut. No constraints are applied to the 
+ p^- and 7r ± functions unless the fits fail. To esti- 
mate the systematic uncertainty of the purity, the mean 
and width of the electron function are allowed to vary 
up to one, two, three and four standard deviations from 
their central values. For each of the four constraints, 
we calculate one value of the purity. The final purity 
is taken as the mean and the systematic uncertainty is 
taken as the largest difference between the mean and the 
four values from the four constraints. To estimate the 
statistical uncertainty of the purity, we rely on a simple 
Monte-Carlo simulation. We first obtain a large sample 
of altered overall na e spectra by randomly shifting each 
data point in the original spectrum in Fig. [5] according to 
a Gaussian distribution with the mean and width set to 
be equal to the central value and the uncertainty of the 
original data point, respectively. We then obtain the pu- 
rity distribution through calculating the purity from each 
of these altered spectra following the same procedure as 
discussed above. In the end, we fit the distribution with 
a Gaussian function and take its width as the statisti- 
cal uncertainty. The total uncertainty of the purity is 
obtained as the quadratic sum of the statistical and sys- 
tematic uncertainties. 

We follow the same procedure in the Run2008 and 
the Run2005 analysis except that the overall na e dis- 
tribution in the Run2008 analysis is the combined re- 
sult from the datasets of all three high-tower triggers as 
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FIG. 7: (a) Derived pr spectrum for inclusive photons (solid 
line) and the uncertainty represented by the region between 
the spectra of 7r° and rj decay photons (dot-dashed line) and 
inclusive photon with doubled direct photon yield (dotted 
line) as well as their ratio to the inclusive photon as shown in 
(b). 
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FIG. 8: (Color online) Photonic electron reconstruction effi- 
ciency as a function of pr for (a) 7 conversion (open circles), 
7T° (closed triangles) and r/ (open triangles) Dalitz decay for 
the Run2008 analysis, (b) combination of 7 conversion, it 
and r\ Dalitz decay for the Run2008 analysis and (c) 7 con- 
version for the Run2005 analysis. The solid line is a fit and the 
dashed lines represent the uncertainty. See text for details. 



described Sec. Ill CI while in the Run2005 analysis, the 
purity are calculated separately for the two high-tower 
triggers. Figure |B] shows the purity as a function of elec- 
tron p T for the Run2008 (a) and the Run2005 (b) data. 
Tighter electron identification cuts and much higher pho- 
tonic electron yield lead to much higher purity for the 
Run2005 inclusive electron sample. 



D. Photonic Electron Reconstruction Efficiency 

Since photon conversions, it and r\ meson Dalitz de- 
cays are the dominant sources of photonic electrons, they 
are the components that we used to calculate e p ho, the 
photonic electron reconstruction efficiency, in the anal- 
ysis of the Run2008 data. The e p h for each individual 
component is calculated separately to account for its pos- 
sible dependence on the decay kinematics of the parent 
particles. The final e p ho is obtained by combining results 
from all components according to their relative contribu- 
tion to the photonic electron yield. 

The determination of e p ho is done through reconstruct- 
ing electrons from simulated 7 conversion or Dalitz decay 
of 7T° and 77 with uniform pt distributions that are em- 
bedded into high-tower trigger events. These events are 
then fully reconstructed using the same software chain as 
used for data analysis. To account for the efficiency de- 
pendence on the parent particle pt, we use a fit function 
to the measured n° spectrum, the derived r\ and inclu- 
sive photon pt spectra as weights. The fit function to the 
measured 7r° spectrum is provided by the PHENIX ex- 
periment in Ref . [32[ . The r\ spectrum is derived from the 
7T° measurement assuming tot scaling, i.e. replacing px 



with y p\ + m 2 h — m^o while keeping the function form 

unchanged. Figure[7](a) shows the derived inclusive 7 px 
spectrum (solid line), and an estimate of its uncertainty 
represented by the region between the dotted and dot- 
dashed lines. Figurc[7](b) shows the uncertainty in linear 
scale. The inclusive 7 spectrum is obtained by adding the 
direct 7 yield to the ir° and r\ decay 7 yield calculated 
using PYTHIA. The direct 7 yield is obtained from the 
fit function to the direct 7 measurement provided by the 
PHENIX experiment in Ref. [32|. The dot-dashed line 
represents the 7 spectrum from 7r° and 77 decay alone. 
The dotted line is obtained by doubling the direct 7 com- 
ponent in the inclusive photon spectrum. By comparing 
the ratio of the derived inclusive 7 yield to that of 7r° 
and r\ decay photon with the double ratio measurement 
in Au+Au most peripheral collisions [33| . we found the 
uncertainty covers the possible variations of the inclusive 
photon yield. 

STAR simulations for 7 conversion and Dalitz decay 
are based on GEANT3 [H[ which incorrectly treats 
Dalitz decays as simple 3-body decays in phase space. 
We therefore modified the GEANT decay routines us- 
ing the correct Kroll-Wada decay formalism (H)]. Their 
kinematics is strongly modified by the dynamic electro- 
magnetic structure arising at the vertex of the transition 
which is formally described by a form factor. We included 
the most recent form factors using a linear approximation 
for the 7T° Dalitz decay [36j . and a pole approximation 
for the decays of rj |37| . 

Figure|8](a) shows the photonic electron reconstruction 
efficiency as a function of electron px for 7 conversion, 
7T° and f] Dalitz decay electrons, which turn out to be 
very similar because of the similar decay kinematics. The 
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increase towards larger electron px is due to the higher 
probability of reconstructing both electrons from high px 
(virtual) photons. The uncertainties shown in the plots 
are dominated by the statistics of the simulated events. 
The effect due to the variation of the inclusive photon 
spectrum shape is found to be negligible for this analy- 
sis. Figure |8] (b) shows the combined photonic electron 
reconstruction efficiency for the Run2008 analysis, which 
is calculated as 

tpho {pt ) = Am „_^ , Zl^lls i ATV^_S ' e 7 (Pt ) + 

^o(pt) 

where and iVJ are respectively the yield of elec- 

trons from photon conversion, 7r° and r\ Dalitz decay; 
e 7 , e n o and e v are the corresponding photonic electron 
reconstruction efficiencies. Based on Table HI approxi- 
mately 36% of the photonic electrons are from ir° Dalitz 
decay and about 10% are from r/ Dalitz decay. Their 
variations have negligible effect on the results since e 7 , 
e^o and are almost identical. The solid line is a second 
order polynomial ht to the data. The systematic uncer- 
tainty is represented by the region between the dotted 
lines, which are second order polynomial fits after mov- 
ing all the data points simultaneously up and down by 
one standard deviation. 

For the Run2005 analysis, the dominant source of pho- 
tonic electrons is conversion in the silicon detectors. We 
therefore neglect contributions from Dalitz decays while 
following the same procedure as for the Run2008 analysis 
to calculate t p ho- Figure [8] (c) shows e p ho as a function 
of pt for 7 conversion for the Run2005 analysis. The 
solid line is a fit to the spectrum with a second order 
polynomial function and the region between dashed lines 
represents the uncertainty estimated in the same way as 
for the Run2008 analysis. The inclusion of the Dalitz 
decays is estimated to reduce the e p ho by less than 0.5% 
which is well within the systematic error. The uncer- 
tainty because of the inaccurate material distribution in 
the simulation as mentioned in Sec. Ill Bl is negligible since 
the majority of the material, dominated by our silicon 
detectors, is within a distance of 30 cm from beam pipe 
and the variation of e p ho of photonic electrons produced 
within this region is small. 



E. Trigger Efficiency 

The trigger efficiency is the ratio of the electron yield 
from high-tower trigger events to that from minimum- 
bias trigger events after normalizing the two according to 
the integrated luminosity. To have a good understand- 
ing of trigger efficiency, one needs enough minimum- 
bias events for the baseline reference. However, for the 
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FIG. 9: (Color online) (a) add) distribution for high-tower 
trigger events, (b) \ 2 as a function of the add) cut. (c) Raw 
inclusive electron px spectrum from VPD trigger in Run2008 
d+Au collisions before (open squares) and after applying the 
add) > 193 cut (closed circles), (d) add) distribution for data 
(closed circles) and simulation (solid line) at px = 4.0 — 5.0 
GeV/c. See text for details. 



Run2008 p + p data, the number of minimum-bias events 
is too small to be used for this purpose. Fortunately the 
Run2008 d+Au data were taken using the same sets of 
high-tower triggers as used for p + p run. Since the two 
data sets were taken serially, the high-tower trigger ef- 
ficiency is expected to be the same. During the rf+Au 
run, many events also were taken using the VPD trigger, 
which is essentially a less efficient minimum-bias trigger 
that can serve as the baseline reference for trigger effi- 
ciency analysis. As a cross check, we also evaluate the 
trigger efficiency through the Run2008 p + p simulation. 

From the VPD trigger events, we first regenerate a 
high-tower trigger px spectrum by requiring add) of 
BEMC clusters to be larger than the threshold. The 
add) is the offline ADC value of a BEMC cluster's most 
energetic tower which is one of the high-towers responsi- 
ble for firing a high-tower trigger. Figure^ (a) shows the 
add) distribution of photonic electrons from high-tower 
trigger events. The sharp cut-off around a value of 200 
is the offline ADC value of the trigger threshold setting. 
The smaller peak below the trigger threshold is due to 
electrons which happen to be in events triggered by some- 
thing else other than the electrons. By requiring add) to 
be larger than the threshold, we reject these electrons 
which did not trigger the event since the uncertainty of 
their yield is affected by many sources and is therefore 
hard to be evaluated reliably. When the threshold is cor- 
rectly chosen, the regenerated spectrum shape should be 
very similar to that of the actual high-tower trigger. We 
therefore quantitatively determine the trigger threshold 
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FIG. 10: (Color online) pr dependence of high-tower trigger 
efficiency from data (closed circles), simulation (closed trian- 
gles) and combined results (solid line) for Run2008 analysis. 
The dashed lines represent the uncertainty. See text for de- 
tails. 
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FIG. 11: (Color online) Raw inclusive electron pr spectrum 
for minimum-bias (closed triangles) and two high-tower trig- 
gers, i.e. HT1 (closed circles) and HT2 (open squares) to- 
gether with a power-law fit (solid line) and fit uncertainty 
(dashed line) for Run2005 analysis. 



as the add) cut which minimizes the 

2 _ V Ni(VPD + adcO) 2 2 
X ^ [ Ni(HT) ' ' a ' 1 

where Ni(VPD + add)) is the regenerated high-tower 
trigger electron yield from VPD events in the i th px bin, 
Ni(HT) is the electron yield at the same px bin from 
the actual high-tower trigger events, and <7j is the uncer- 
tainty of N q (VPD + add)). Figure M 0) shows the \ 2 
as a function of the add) cut; the threshold is taken as 
193. Figure [9] (c) shows the px spectrum of raw inclu- 
sive electrons from the VPD trigger (open squares) and 
the regenerated high-tower spectrum (closed circles) af- 
ter applying the add) > 193 cut used to calculate the 
trigger efficiency. 

To estimate trigger efficiency through simulation, we 
tune the simulated single electron add) spectrum in each 
individual px bin to agree with the data in the region 
above the threshold. The data spectra are obtained from 
the unlike-minus-like pairs, i.e. pure photonic electrons. 
As a demonstration of the comparison, Fig. [9] (d) shows 
the spectra from data (closed circles) and simulation 
(solid line) at 4.0 GeV/c < p T < 5.0 GeV/c. The ef- 
ficiency is defined as the fraction of the simulated add) 
spectrum integral above the trigger threshold. 

In the Run2008 analysis the raw px spectrum of non- 
photonic electrons is obtained by combining the datasets 
of all three high-tower triggers. Since the shape of the 
combined spectrum is the same as that of the high-tower 
trigger with the lowest threshold, we only need to esti- 
mate the trigger efficiency of this lowest threshold trigger. 
Figure [TO] shows the trigger efficiency as a function of px 
that is calculated using d+Au VPD events (closed circles) 
and simulated events (closed triangles) in the Run2008 
analysis. At px > 3.5 GeV/c, they agree with each other 
reasonably well. At lower px, the simulated results are 



not reliable because the numerator in the efficiency cal- 
culation is only from a tail of the spectrum and a small 
mismatch between simulation and data can have a large 
impact on the results. On the other hand, the results 
from VPD events suffer from low statistics at high px- 
The final efficiency is therefore taken to be the combina- 
tion of the two, i.e., at pt < 3.5 GeV/c, the efficiency 
is equal to that from VPD events assigning a systematic 
uncertainty identical to the statistical uncertainty of the 
data point, while at high px, the efficiency is equal to the 
simulated results, and the systematic error are from the 
tuning uncertainty. 

In the Run2005 analysis, the efficiencies of the two 
high-tower triggers are estimated separately. While there 
arc more minimum-bias events for Run2005 than for 
Run2008, the statistics arc poor at px > 2.0 GeV/c. 
We thus rely on a fit to the spectrum, which consists 
of minimum-bias events at low px and high-tower trig- 
ger events at high px where the trigger is expected to 
be fully efficient, as the baseline reference for the trigger 
efficiency evaluation. Figure [TT] shows the raw inclusive 
electron px spectrum from minimum-bias, HT1 and HT2 
events. The fit uses a power-law function A(l+px/B)~ n . 
The regions where we expect HT1 and HT2 trigger to be 
fully efficient are above 4.5 and 6.0 GeV/c respectively. 
The dashed line shows the fit uncertainty, obtained from 
many fit trials. In a single fit trial, each data point is 
randomized with a Gaussian random number, with the 
mean to be the central value and the rms to be the statis- 
tical uncertainty of the data point. Additional systematic 
uncertainty coming from fits using different functions is 
included in the cross section calculation and is not dis- 
played in the figure. Figure [12] shows the efficiency of 
HT1 (a) and HT2 (b) triggers, defined as the ratio of 
the raw HT1 or HT2 inclusive electron spectrum to the 
baseline fit function. We used error functions to param- 
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FIG. 12: (Color online) pr dependence of trigger efficiency 
for the two high-tower triggers, i.e. HT1 (a) and HT2 (b) 
including result from data (closed circles) for Run2005, an 
error function fit (solid lines) and the fit uncertainty (dot- 
dashed lines). 



FIG. 13: (Color online) (a) Distribution of the calculated 
ubbc before (solid lines) and after (closed circles) removing 
events at the beginning and the end of Run2008. (b) Variation 
of the calculated obbc a.s a function of run number. The runs 
outside the region between the two dashed lines are rejected. 



eterize both efficiencies. The dashed lines represent the 
uncertainty, obtained in the same way as for the Fig. 1111 
fits. 

High-tower trigger efficiency for photonic and non- 
photonic electrons can be different. Unlike non-photonic 
electrons, a photonic electron always has a partner. In 
case both share the same BEMC tower, the deposited 
energy will be higher than that for an isolated electron 
and will lead to a higher efficiency. The effect can be 
quantified by comparing the ratio of the isolated electron 
yield to the yield of electrons with partners in minimum- 
bias events to the same ratio in high-tower trigger events. 
We found that the difference is negligible at pt > 2.5 
GeV/c, while the trigger efficiency for photonic electron 
is 20 — 30% higher than for non-photonic electron in the 
lower pt region. 



minosity in Run2008. The minor one centered at a higher 
value comes from events taken at the beginning and the 
end of Run2008 represented by the regions beyond the 
two dashed lines in Fig. [13] (b) showing the calculated 
BBC cross section as a function of run number. After 
removing these runs taken at the beginning and the end 
of Run2008, the minor peak disappeared and the perfor- 
mance of the monitor appeared to be very stable. 

In the data analysis, we also reject those with ctbbc < 
20 mb or ubbc > 30 mb. We fit the <?bbc distribution 
with a Gaussian function and assign the width of the 
function (2.3%) as the systematic uncertainty of the a mon 
with respect to BBC cross section. 

The integrated luminosity sampled by the high-tower 
triggers are ^2.6 pb _1 and ~2.8 pb _1 for Run2008 and 
Run2005, respectively. 



F. Stability of the Luminosity Monitor 

The BBC trigger was used to monitor the integrated 
luminosity for Run2005. During Run2008, because of the 
large beam background firing the BBC trigger, a high 
threshold high-tower trigger was used as the luminosity 
monitor. To quantify the stability of the monitor with 
respect to BBC, we calculate the BBC cross section as a 
function of run number using obbc = N m i n i>i as /£, where 
C = N mon /a mon , a mon is the monitor cross section which 
is estimated to be 1.49 /jb using low luminosity runs, 
Nminbias and N mon are respectively the number of events 
from the BBC trigger and the monitor after correcting 
for prescaling during data acquisition. Here a run refers 
to a block of short term (^30 minutes) data taking. Fig- 
ure [13] (a) shows the distribution of the calculated (Jbbc- 
There are two peaks in the figure. The dominant one 
centered around 25 mb contains most of the recorded lu- 



G. Contribution from Vector Mesons 

The main background sources of electrons that do not 
originate from photon conversion and Dalitz decay are 
electromagnetic decays of heavy {J/ip, T) and light vec- 
tor mesons (p, u> and <f>) as well as those from Drell-Yan 
process. 

The electrons from J decay contribute noticeably to 
the observed non-photonic electron signal as pointed out 
in Ref. [38[. In order to estimate the contribution from 
J/ip —> e + e~ to the non-photonic electron yield, we com- 
bine the measured differential J/ip cross-sections from 
PHENIX fH and STAR [H. For each data point we 
add the statistical and systematic uncertainties, except 
the global uncertainties, in quadrature. Figure [14] (a) 
shows the measured J/ip differential cross section from 
the two experiments. While the PHENIX measurement 
dominates the low to medium-p^ region, the STAR mea- 
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surcmcnt dominates the high-p^ region. The combined 
spectrum is fit using a power-law function of the form 
Ed 3 a/d 3 p\ y = = A(exp(apT — bp T 2 ) +Pr/po)~ n where 

a = 0.32±0.04 GeV" 1 ^ b 
po = 2.59±0.21 GeV/c and n = 
8.44±0.61 are fit parameters. The x 2 / NDF 01 the fit is 
27.8/25. To obtain the uncertainty of the fit, the global 
uncertainties of the STAR and the PHENIX (10% jU) 
measurements are assumed to be uncorrelated. We move 
the PHENIX data up by 10% and repeat the fit to obtain 
the band of 68% confidence intervals. The upper edge of 
the band is treated as the upper bound of the fit. Fol- 
lowing the same procedure except moving the PHENIX 
data down by 10%, we obtain the lower bound of the fit 
as the lower edge of the band. Furthermore, since we are 
considering a rather large pr range (pr < 14.0 GcV/c ), 
we cannot assume that the pt and rapidity distributions 
factorize. We use PYTHIA to generate dN/dy(pT) and 
implement a Monte-Carlo program using the above func- 
tions as probability density functions to generate J/ij) 
and decay them into e + e~ assuming the J/tp to be un- 
polarized. The decay electrons are filtered through the 
same detector acceptance as used for the non-photonic 
electrons. The band in Fig. [14] (b) shows the invariant 
cross section of J/ip decay electrons as a function of the 
electron pr- The uncertainty of the derived yield comes 
from the uncertainty of the fit to the J /if; spectra and is 
represented by the band which is also shown in Fig. [14] 
(c) in linear scale. 

The invariant cross section of electrons from T de- 
cays (T — ¥ e + e~), represented by the dot-dashed line in 
Fig. HH (b), is calculated in a similar fashion as that for 
the J/ij) except that the input T spectrum is from a Next- 
to-Leading Order pQCD calculation in the color evapora- 
tion model (CEM) [42(. We have to rely on model calcu- 
lations since so far no invariant px spectrum in our energy 
range has been measured. However, in a recent measure- 
ment STAR reported the overall production cross section 
for the sum of all three T(1S+2S+3S) states in p+p colli- 
sions at yfs = 200 GcV tobcBx dajdy = 114± 38^4 pb, 
which is consistent with the CEM prediction [43[ . Adding 
the statistical and systematic uncertainty in quadrature, 
the total relative uncertainty of this measurement is 
^39%, which is the value we assigned as the total un- 
certainty of the T feed-down contribution to the non- 
photonic electrons at all px- 

The contribution to the non-photonic electron yield 
from the light vector mesons is estimated using PYTHIA, 
assuming the meson spectra follow tut scaling. We gener- 
ate a sample of decay electrons using light vector mesons 
with flat spectra in pt as input. To derive the differential 
cross section of the electrons, we keep only those electrons 
within the same detector acceptance as that for the non- 
photonic electrons and weight them with the spectra of 
p, lu and <f>. The meson spectra are obtained by replacing 

the pt with 



P T 



ml 



m r , in the same fit function as 
for the 7T° measurement (see Fig. [7]). Here is the mass 
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FIG. 14: (Color online) (a) The J/ij) invariant cross section 
measurement from STAR (closed circles) and PHENIX (open 
triangles), together with the fits using A(exp(apr — bp^) + 
Pt/po) " (solid line), (b) Invariant cross section of the elec- 
tron from decays of J/ij) (band), T (dot-dashed line), Drell- 
Yan (dotted line) and light vector mesons (solid line). The 
uncertainty of the J/i/j feed-down is represented by the band 
shown in (c) in linear scale. 



of the vector meson 

tO 7T 



The relative yields of the mesons 
17] are also taken into account during this process. 
We include the decay channels <\> — > e + e~, <\> — > rje + e~, 
uj — > e + e~ , ui — > 7r°e + e~ and p — > e + e~ in the calcu- 
lation. The derived electron differential cross section is 
represented by the solid line in Fig. [T?] (b). We assign a 
50% systematic uncertainty to cover the uncertainty of 
the 7T° measurement and the meson to pion ratios. 

The contribution to the non-photonic electron yield 
from the Drell-Yan processes is represented by the dotted 
line in Fig. [14] (b) and is estimated from a Leading-Order 
pQCD calculation using the CTEQ6M parton distribu- 
tion function with a K-factor of 1.5 applied and without 
a cut on the electron pair mass [44) ]. No uncertainty is 
assigned to this estimate. 



IV. RESULTS 

A. Non-photonic Electron Invariant Cross Section 

The invariant cross section for non-photonic electron 
production is calculated according to 



E- 



4* a 
dp 3 
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N„ 
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£ 2np T Ap T Ay 

tree ttrig ^eid ^BBC 



where N npe is the non-photonic electron raw yield with 
the V z cuts, e rec is the product of the single electron re- 
construction efficiency and the correction factor for mo- 
mentum resolution and finite spectrum bin width, ttrig is 
the high-tower trigger efficiency, e e id is the electron iden- 
tification efficiency, C is the integrated luminosity with 
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TABLE II: Sources of systematic uncertainty for the non- 
photonic electron invariant yield in p+p collisions. Type A are 
point to point uncertainties. Type B are scaling uncertainties 
which move data points in the same direction. Type C are 
the scaling uncertainties that are common to both Run2008 
and Run2005. The range in each individual source covers the 
variation of the systematic uncertainty as a function of pr- 



source 


Run2008 


Run2005 




5.0-48.1 % (A) 


8.5-38.0 % (A) 




6.5-25.2 % (A) 


0.7-2.0 % (A) 


e trg 


1.8-10.0 % (A) 


0.3-16 % (A) 




5.4 % (B) 




f-rec 


2.3-33.3 % (A) 


1.0-3.5 % (A) 




15.7 % (B) 


11.0 % (B) 


C-£BBC 


2.3 % (B) 






8.1 % (C) 


8.1 % (C) 



the V z cuts and e B BC = 0.866 ± 0.08 is the BBC trigger 
efficiency. The systematic uncertainties of all these quan- 
tities are listed in Tabic [TTJ The relative uncertainty of 
£'£bbc m maximum r ang e is 14% including uncertainties 
in tracking efficiency |24j. Assuming a flat distribution 
within the range, we estimate the C-cbbc uncertainty to 
be 8.1% in one standard deviation. The uncertainty of 
N npe is the quadratic sum of the uncertainties from the 
estimation of e p ho, purity and the light vector meson con- 
tribution. The uncertainty of e rec is the quadratic sum 
of the uncertainties from correcting the track momentum 
resolution, the finite spectrum bin width as well as the es- 
timation of single electron reconstruction efficiency. The 
range of the uncertainty for each individual quantity cov- 
ers the variation of the uncertainty as a function of pt ■ 
In order to compare with the result in Ref. [HI, Qj} , we do 
not subtract the J/ip> T an d Drell-Yan contribution from 
the non-photonic electron invariant cross section shown 
in Fig. COD and Fig. [16] 

Figure [15] (a) shows the ratio of non-photonic to pho- 
tonic electron yield as a function of px in p + p collisions 
in Run2008 (closed circles) and Run2005 (open trian- 
gles). The ratio for Run2008 is much larger because there 
was much less material in front of the TPC for Run2008. 
Figure [15] (b) shows the non-photonic electron invariant 

cross section ( e t e ) as a function of pr in p + p colli- 
sions from the Run2008 analysis (closed circles) and the 
Run2005 analysis (open triangles). Despite the large dif- 
ference in photonic background, the two measurements 
are in good agreement. 

Figure [16] (a) shows the non-photonic electron ( e + e ) 
invariant cross section obtained by combining the 
Run2008 and the Run2005 results using the "Best Lin- 
ear Unbiased Estimate [45| . The corrected result of our 
early published measurement using year 2003 data [l6| is 
shown in the plot as well. The published result exceeded 
pQCD predictions from FONLL calculations by about a 
factor of four. We, however, uncovered a mistake in the 



corresponding analysis in calculating e p ho- The details 
are described in the erratum [l6j]. To see more clearly 
the comparison, Fig. [16] (b) shows the ratio of each indi- 
vidual measurement, including PHENIX results, to the 
FONLL calculation. One can see that all measurements 
at RHIC on non-photonic electron production in p+p col- 
lisions are now consistent with each other. The corrected 
run 2003 data points have large uncertainties because 
of the small integrated luminosity (^100 nb" 1 ) in that 
run. FONLL is able to describe the RHIC measurements 
within its theoretical uncertainties. 



B. Invariant Cross Section of Electrons from 
Bottom and Charm Meson Decays 

Electrons from bottom and charm meson decays are 
the two dominant components of the non-photonic elec- 
trons. Mostly due to the decay kinematics, the azimuthal 
correlations between the daughter electron and daughter 
hadron are different for bottom meson decays and charm 
meson decays. A study of these azimuthal correlations 
has been carried out on STAR data and is compared 
with a PYTHIA simulation to obtain the ratio of the 
bottom electron yield to the heavy flavor decay electron 
yield {e B /{e B + e D )) [Til, where PYTHIA was tuned to 
reproduce STAR measurements of D mesons px spectra 
[46j |. Using the measured e_B/(ee + erj), together with 
the measured non-photonic electron cross section with 
the electrons from J/i/j, T decay and Drell-Yan processes 
subtracted, we are able to disentangle these two compo- 
nents. 




2468 2468 10 



P T (GeV/c) p T (GeV/c) 

FIG. 15: (Color online) (a) Ratio of non-photonic to pho- 
tonic electron yield from the Run2008 (closed circles) and the 
Run2005 (open triangles) analyses, (b) Invariant cross sec- 
tion for non-photonic electron production ( e + e ) in p + p 
collisions from the Run2008 (closed circles) and the Run2005 
(open triangles) analyses. The error bars and the boxes repre- 
sent statistical and systematic uncertainty, respectively. The 
solid line is FONLL calculation and the dashed lines are the 
FONLL uncertainties 0. 
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FIG. 16: (Color online) (a) Invariant cross section of non- 
photonic electron production (- 



-) in p + p collisions 
from this analysis (closed circles) after combining results 
from Run2005 and Run2008. The published STAR re- 
sult [l6l [ (closed triangles) is also shown, (b) Ratio of data 
over FONLL 0] from all measurements at RHIC including 
PHENIX results [l7| (open triangles). 



The bottom electron cross section is calculated as 
es/(e_B + en) times the non-photonic electron cross sec- 
tion with the contribution from J/ip, T decay and Drell- 
Yan processes subtracted. The same procedure applies 
to the charm electrons except that (1 — €3/(^3 + ejj)) 
is used instead. The specific location of px where the 
£-b/{sb + gd) is measured, is different from that of the 
non-photonic electrons. To accommodate the difference, 
we calculate es/(es + en) in any given px in non- 
photonic electron measurements through a linear inter- 
polation of the actual es/(es + ec) measurements. As 
an estimation of the systematic uncertainty of the inter- 
polated value, we also repeat the same procedure using 
the curve predicted by FONLL. Figure [T71 shows the in- 
variant cross section of electrons (— ^—) from bottom 
(upper-left) and charm (upper-right) mesons as a func- 
tion of pt and the corresponding FONLL predictions, 
along with the ratio of each measurement to the FONLL 
calculations (lower panels). The statistical uncertainty of 
each data point is obtained by adding the relative statis- 



tical uncertainties of the corresponding data points in the 
non-photonic electron and the e# / '(es+eu) measurement 
in quadrature. The systematic uncertainties are treated 
similarly, except that the uncertainties from the interpo- 
lation process are also included. The measured bottom 
electrons are consistent with the central value of FONLL 
calculation and the charm electrons arc in between the 
central value and upper limit of the FONLL calculation, 
the uncertainties of which are from the variation of heavy 
quark masses and scales. From the measured spectrum, 
we determine the integrated cross section of electrons 
at 3 GeV/c < p T < 10 GcV/c from bottom 
and charm meson decays to be, respectively, 



da 



(B-te)4-(B^fl^e) 

dy e 



da j 



dy e 



ly«=o 



4.0 ± 0.5(stat.) ± l.l(syst.)nb 



6.2 ± 0.7(stat.) ± 1.5(syst.)nb, 



where y e is the electron rapidity. The 8.1% global scale 
uncertainty from the BBC cross section is included in the 
total systematic uncertainty. 

Relying on theoretical model predictions to extrapo- 
late the measured results to the phase space beyond the 
reach of the experiment, one can estimate the total cross 
section for charm or bottom quark production. We per- 
form a PYTHIA calculation with the same parameters 
as in Ref. [lj| . After normalizing the pr spectrum to our 
high-pT measurements and extrapolating the results to 
the full kinematic phase space, we obtain a total bottom 
production cross section of 1.34 /A>. However, with the 
PYTHIA calculation using the same parameters except 
MSEL— 5, i.e. bottom production processes instead of 
minimum-bias processes as in the former calculation, we 
obtain a value of 1.83 /jb. The PYTHIA authors rec- 
ommend the minimum-bias processes [3l| . This large 
variation between the extracted total bottom produc- 
tion cross sections comes mostly from the large differ- 
ence in the shape of the bottom electron spectrum in 
the two PYTHIA calculations with MSEL=1 and with 
MSEL=5. Since both calculations are normalized to 
the measured data, the difference in the shape shows 
up at p T < 3 GeV/c. The fact that the PYTHIA cal- 
culation with MSEL=5 only includes leading order di- 
agrams of bottom production causes the difference be- 
tween the PYTHIA calculations. Measurements in the 
low pt region are therefore important for the understand- 
ing of bottom quark production at RHIC. Both values are 
consistent with the FONLL @ prediction, imt^^h, 
within its uncertainty. 



V. CONCLUSIONS 

STAR measurements of high pt non-photonic electron 
production in p + p collisions at ^/i= 200 GeV using 
data from Run2005 and Run2008 agree with each other 



16 




FIG. 17: (Color online) Invariant cross section of electrons ( e ) from bottom (upper-left) and charm meson (upper-right) 
decay, together with the ratio of the corresponding measurements to the FONLL predictions for bottom (lower-left) and charm 
electrons (lower-right). The solid circles are experimental measurements. The error bars and the boxes are respectively the 
statistical and systematic uncertainties. The solid and dotted curves are the FONLL predictions and their uncertainties. The 
dashed and dot-dashed curves are the FONLL prediction for B— >D— >e, i.e. electrons from the decays of D mesons which in 
turn come from B meson decays. 



despite the large difference in background. This mea- 
surement and PHENIX measurement are consistent with 
each other within the quoted uncertainties. After correct- 
ing a mistake in the photonic electron reconstruction effi- 
ciency, the published STAR result using year 2003 data is 
consistent with our present measurements. We are able 
to disentangle the electrons from bottom and charm me- 
son decays in the non-photonic electron spectrum using 
the measured ratio of es/(es + gd) and the measured 
non-photonic cross section. The integrated bottom and 
charm electron cross sections ( e j, 6 ) at 3 GeV/c < pr < 
10 GeV/c are determined separately as 



da 



(B-te)+(B-tD-te) 



dap- 

dy e 



4.0 ± 0.5(stat.) ± l.l(syst.)nb 



,o = 6.2 ± 0.7(stat.) ± 1.5(syst.)nb. 



FONLL can describe these measurements within its 
theoretical uncertainties. Future measurements on low- 
Pt electrons from bottom meson decay are important to 



overcome the large uncertainties of the derived total bot- 
tom quark production cross section that originate mostly 
from the large variations of theoretical model prediction 
in the low-pr region. 
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